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(^■^ ' Theoretical foundations and applications of the generalized spin-fermion ( sp-d ) exchange lattice 

I model to various magnetic systems, e.g. rare-earth metals and compounds and magnetic semicon- 

, ductors are discussed. The capabilities of the model to describe spin quasi-particle spectra are 

investigated. The main emphasis is made on the dynamical behavior of two interacting subsystems, 
the localized spins and spin density of itinerant carriers. A nonperturbative many-body approach, 
the Irreducible Green Functions (IGF) method, is used to describe the quasi-particle dynamics. 
Scattering states are investigated and three branches of magnetic excitations are calculated in the 
regime characteristic of a magnetic semiconductor. For a simplified version of the model ( Kondo 
lattice model ) we study the spectra of quasi-particle excitations with special attention given to 
diluted magnetic semiconductors in simple approximation, to demonstrate the possibilities of the 
IGF approach. For this, to include the effects of disorder, a modified mean fields are determined 
^ ] self-consistently. The approach permits to investigate and clarify the role of various interactions 

I . and disorder effects in unified and coherent fashion. 

^ ■ PACS numbers: Tl.lOFd, Tl.lOLi, 75.10.-b, 75.10Jm, 75.10Lp, 75.30.-m, 75.30.Ds, 75.50.Pp 
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I. INTRODUCTION 



B 

^ ' Existence and properties of localized and itinerant magnetism in metals, oxides and alloys and their interplay is an 
Q interesting but not yet fully understood problem of quantum theory of magnetismiiSi^. Behavior and the true nature 
Q ■ of the electronic and spin states, and their quasi-particle dynamics are of central importance to the understanding of 
physics of correlated systems such as magnetism and Mott-Hubbard metal- insulator transition in metals and oxides, 
magnetism and heavy fermions (HF) in rare-earths compounds, and anomalous transport properties in perovskite 
^ manganites. This class of systems is characterized by the complex, many-branch spectra of elementary excitations. 

■ Moreover, the correlations effects ( competition and interplay of Coulomb correlation, direct or indirect exchange, 
' sp-d hybridization, electron-phonon interaction, disorder, etc.) are essential^. These materials are systems of great 

interest both intrinsically and as a possible source of understanding the magnetism of matter generally. Beginning 
25 with Zener^'^'*'^'^°, De Gennesi^, Ruderman and Kittel^"'^, and Doniach^^, various formulations of spin-fermion model 
for the interacting spin and charge subsystems have been studied. There has been considerable interest in identifying 
the microscopic origin of quasi-particle states in such systems and a few model approaches have been proposed. Many 
magnetic and electronic properties of rare-earth metals and compounds-^ and magnetic semiconductorsMii^ and re- 

■ lated materials may be interpreted reasonably in terms of combined spin-fermion models ( SFM ) which include the 
' interacting spin and charge subsystems2i2iiiiSiiLi&. This approach permits one to describe significant and interesting 

I physics, e.g., the bound states and magnetic polaronsi^fSS, anomalous transport properties, etc. 
■ O ' The problem of adequate physical description within various types of generalized spin-fermion model has intensively 
C , been studied during the last decades, especially in the context of magnetic and transport properties of rare-earth and 
^ ' transition metals and their compoundsi;^ and magnetic semiconductorsiSii2iS£. 

More recent efforts have been directed to the study of the properties of magnetioi^*i^ and diluted magnetic 
semiconductor9^ii22iS2iSli2Si2Si21. This field is very active and there are many aspects to the problem. A lot of materials 
were synthesized and tested22i22i22i2i. The new materials design approach to fabrication of new functional diluted 
magnetic semiconductors (DMS) resulted in producing a variety of compounds . Diluted magnetic semiconductors 
are semiconducting alloys whose lattice contains magnetic atoms as randomly distributed substitutional impurities 
such as Mn-doped InAs or GaAs ( general formula A{^J^MnxB^ ). A fraction of A sublattice which is substituted 
at random by Mn changes the carrier density from the low doping ( insulating regime) to the large doping (metallic 
regime ) . The presence of the spin degree of freedom in DMS may lead to a new semiconductor spin electronics which 
will combine the advantages of the semiconducting devices with the new features due to the possibilities of controlling 
the magnetic state. 

However, the coexistence of ferromagnetism and semiconducting properties in these compounds require a suitable 
theoretical model which would describe well both the magnetic cooperative behavior and semiconducting properties 
as well as a rich field of interplay between them. The majority of theoretical papers on DMS studied its properties 
mainly within the mean field approximation and continuous media terms. In such a picture the disorder effects, which 
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play an essential role'^2i2iMi25i^, can be taken into account, as a rule, roughly only. Moreover, there are different 
opinion on the intrinsic origin and the nature of disorder in DMS^£*2Si^. Recently, there were made a lot of efforts to 
go beyond the simplest level of approximation, the virtual crystal approximation (VGA) and many effective schemes 
for better threatment of disorder effects were elaboratedi22i^ii2Sii2iiiii2iM44 (for a detailed review see Refsi22i2^). Thus 
many experimental and theoretical investigations call for a better understanding of the relevant physics and nature of 
solutions ( especially magnetic ) in terms of the lattice spin-fermion model^'^^'^S'^. In this paper we will concentrate 
on the description of the magnetic excitation spectra and will threat the disorder effects in simplest VGA to emphasize 
the chief purpose of this paper, the need for suitable definition of the relevant generalized mean fields and for internal 
self-consistency in the description of the spin quasiparticle many-body dynamics. 

In the previous papers, we set up the formalism of the method of Irreducible Green Functions (IGF)^. This IGF 
method allows one to describe the quasi-particle spectra with damping for a many-particle system on a lattice with 
complex spectra and a strong correlation in a very general and natural way. This scheme differs from the traditional 
method of decoupling of an infinite chain of equations^ and permits construction of the relevant dynamical solutions 
in a self-consistent way at the level of the Dyson equation without decoupling the chain of equations of motion for 
the GFs. 

In this paper, we apply the IGF formalism to consider the quasi-particle spectra for the lattice spin-fermion model 
consisting of two interacting subsystems. It is the purpose of this paper to explore more fully the notion of General- 
ized Mean Fields (GMF)& which may arise in the system of interacting localized spins ( including elfects of disorder ) 
and lattice fermions to justify and understand the "nature" of the relevant mean fields. Theoretical foundations and 
applications of the generalized spin-fermion ( sp-d ) exchange model to magnetic and diluted magnetic semiconductors 
are discussed in some detail. The capabilities of the model to describe quasi-particle spectra are investigated. The 
key problem of most of this work has remained the formation of spin excitation spectra under various conditions on 
the parameters of the model. The intention is to investigate the quasi-particle spectra and GMF of the magnetic 
semiconductors consisting of two interacting charge and spin subsystems within the lattice spin-fermion model in a 
unified and coherent fashion to analyze the role and influence of the Goulomb correlation and exchange. The chief 
purpose of this paper has been to call attention to the need for internal self-consistency in the description of spin 
quasi-particle dynamics of interacting spin and charge subsystems. 



II. SPIN-FERMION MODEL 



The concept of the sp — d ( or d — f ) model plays an important role in the quantum theory of 
magnetismi*2i2i^iiSii2iiSii&. In this section, we consider the generalized sp ~ d model which describes the localized 
3(i(4/)-spins interacting with s(p)-like conduction (itinerant) electrons ( or holes ) and takes into consideration the 
electron-electron interaction. 
The total Hamiltonian of the model is given by 

H = H,+H,^d + Hd (1) 
The Hamiltonian of band electrons (or holes) is given by 

Hs = J2Yl ^»j"Laj'T + n^^n.^-^ (2) 

ij <7 i{7 

This is the Hubbard model. We adopt the notation 

a^a = N-^/^ aka exp{ikR^) = N'^/^ ^ 4^ exp(-ifci?i) 

k k 

In the case of a pure semiconductor, at low temperatures the conduction electron band is empty and the Goulomb 
term U is therefore not so important. A partial occupation of the band leads to an increase in the role of the Goulomb 
correlation. It is clear that we treat conduction electrons as s-electrons in the Wannier representation. In doped DMS 
the carrier system is the valence band p-holes. 
The band energy of Bloch electrons e(fc) is defined as follows: 

t^J = ^ e(fc) exp[ik{R, - Rj)], 

k 
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where N is the number of lattice sites. For the tight-binding electrons in a cubic lattice we use the standard expression 
for the dispersion 

e(fc) = 2'y^ t{da) cos{kaa) (3) 
, where Oq, denotes the lattice vectors in a simple lattice with the inversion centre. 

The term Hs-d describes the interaction of the total 3d(4f)-spin with the spin density of the itinerant carriers^ 

Hs-d = -2 E ^'^^^^ = -IN-'/^ E E[^->L«fc+«-. + z,Sl^alak+,.] (4) 

i kq (7 

where sign factor z„ is given by 

Zcr = (+or-) for (7 = (t or J,) 

and 

^-a ^ ( ^-q if = + 
-1 \Sl^ if CT = - 

In DMS the local exchange coupling resulted from the p — d hybridization between the Mn d levels and the p valence 
band / ^ V^_^ . For the subsystem of localized spins we have 



^d, — r, E/ Jij^i^'j ~ r, E/ Jq^q^-q (5) 

Here we use the notation 



2^ ' ^ 2 



Sr = N-'/^ E eM^kR^) = N-'/^ E exp(-ifci?.0 

k i 

± 2 ± 2 

k 

This term describes the direct exchange interaction between the localized 3d (4f) magnetic moments at the lattice 
sites i and j. In the DMS system this interaction is rather small. The ferromagnetic interaction between the local 
Mn moments is mediated by the real itinerant carriers in the valence band of the host semiconductor material. The 
carrier polarization produces the RKKY exchange interaction of Mn local moments 

Hrkky — —^"^KijSiSj (6) 

We emphasize that Kij ^ |/^| ^ Vp_j^. To explain this, let us remind that the microscopic model^, which contains 
the basic physics, is the Anderson-Kondo model 

iJ = E E *y"L«J'^ -V^^ia+dja + h.c.) 

ij a ij a 

For the symmetric case U = 2Ed and for ^ y Eq.Q can be mapped onto the Kondo lattice model ( KLM ) 

= E E - E ^^'^'^^ (8) 



IJ CT 



Here I ^ The KLM may be viewed as the low-energy sector of the initial model Eq. ( [T] ) 
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III. OUTLINE OF THE IGF METHOD 



In this section, we discuss the main ideas of the IGF approach that aUows one to describe completely quasi-particle 
spectra with damping in a very natural way. 

We reformulated the two-time GF method^ to the form which is especially adjusted to correlated fermion systems on 
a lattice and systems with complex spectra^i^*i&. A very important concept of the whole method is the Generalized 
Mean Fields (GMFs), as it was formulated in^. These GMFs have a complicated structure for the strongly correlated 
case and complex spectra, and are not reduced to the functional of mean densities of the electrons or spins when one 
calculates excitation spectra at finite temperatures. 
To clarify the foregoing, let us consider a retarded GF of the form^ 

=« A{t),A\t') »= -te{t - t') < [A{t)A\t% >,T]^±l (9) 

As an introduction to the concept of IGFs, let us describe the main ideas of this approach in a symbolic and simplified 
form. To calculate the retarded GF G{t — t'), let us write down the equation of motion for it 

ivGiLu) =< [A, A% > + « [A, H]^ I At »^ (10) 

Here we use the notation << A{t), A^ {t') >> for the time-dependent GF and << A \ A'' >>lu for its Fourier 
transform^S. The notation [A, B]jj refers to commutation and anticommutation depending on the value of 77 = ±. 
The essence of the method is as follows^: 

It is based on the notion of the "IRREDUCIBLE" parts of GFs (or the irreducible parts of the operators, A and A^ , 
out of which the GF is constructed) in terms of which it is possible, without recourse to a truncation of the hierarchy 
of equations for the GFs, to write down the exact Dyson equation and to obtain an exact analytic representation for 
the self-energy operator. By definition, we introduce the irreducible part (ir) of the GF 

« [A, H]_\A^ »=« [A, H]_- zA\A^ » (11) 

The unknown constant z is defined by the condition (or constraint) 

<[[A,H]^^\A%,>=0 (12) 

which is an analogue of the orthogonality condition in the Mori formalism. From the condition H12I) one can find 

_ <[[A,H].,AX> _Ah 

<[AA%> -Mo ^'^^ 

Here Mq and Mi are the zeroth and first order moments of the spectral density. Therefore, the irreducible GFs are 
defined so that they cannot be reduced to the lower-order ones by any kind of decoupling. It is worth noting that 
the term "irreducible" in a group theory means a representation of a symmetry operation that cannot be expressed in 
terms of lower dimensional representations. Irreducible (or connected ) correlation functions are known in statistical 
mechanics. In the diagrammatic approach, the irreducible vertices are defined as graphs that do not contain inner 
parts connected by the G'^-line. With the aid of the definition l|ll() these concepts are expressed in terms of retarded 
and advanced GFs. The procedure extracts all relevant (for the problem under consideration) mean-field contributions 
and puts them into the generalized mean-field GF which is defined here as 

^0 ^ <[AAX > 
(w - z) 

To calculate the IGF « [A, H]^{t), A"' {t') >> in (irU)) . we have to write the equation of motion for it after 

differentiation with respect to the second time variable t' . The condition of orthogonality p2|l removes the inhomo- 
geneous term from this equation and is a very crucial point of the whole approach. If one introduces the irreducible 
part for the right-hand side operator as discussed above for the "left" operator, the equation of motion (fTU)) can be 
exactly rewritten in the following form 

G = G" + G"PG° (15) 

The scattering operator P is given by 

P= (Mo)-i( (-) « [A,HU[AlH]^ »(^'^))(Afo)-' (16) 
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The structure of equation ( I16|) enables us to determine the self-energy operator M, by analogy with the diagram 
technique 

P = M + MG"P (17) 

We use here the notation M for self-energy ( mass operator in the quantum field theory ). From the definition (|17|l it 
follows that the self-energy operator M is defined as a proper (in the diagrammatic language, "connected") part of the 
scattering operator M = (Py. As a result, we obtain the exact Dyson equation for the thermodynamic double-time 
Green functions 

G = G" + G°MG (18) 

The difference between P and Af can be regarded as two different solutions of two integral equations (|15|l and H18|l . 
But from the Dyson equation (|18() only the full GF is seen to be expressed as a formal solution of the form 

G^ [(G°)-i (19) 

Equation (|19|l can be regarded as an alternative form of the Dyson equation p8|l and the definition of M provided 
that the generalized mean-field GF G° is specified. On the contrary , for the scattering operator P, instead of the 
property G^G^^ + G^M = 1, one has the property 

(G°)-^ - G-i = PG°G^i 

Thus, the very functional form of the formal solution H19|) determines the difference between P and M precisely. 
Thus, by introducing irreducible parts of GF (or irreducible parts of the operators, out of which the GF is constructed) 
the equation of motion (|10() for the GF can exactly be ( but using the orthogonality constraint 1^21) transformed into 
the Dyson equation for the double-time thermal GF (|18|l . This result is very remarkable because the traditional form 
of the GF method does not include this point. Notice that all quantities thus considered are exact. Approximations 
can be generated not by truncating the set of coupled equations of motions but by a specific approximation of the 
functional form of the mass operator M within a self-consistent scheme expressing M in terms of initial GF 

M « F[G] 

Different approximations are relevant to different physical situations. 

The projection operator technique has essentially the same philosophy. But with using the constraint (|12|l in our 
approach we emphasize the fundamental and central role of the Dyson equation for the calculation of single-particle 
properties of many-body systems. The problem of reducing the whole hierarchy of equations involving higher-order 
GFs by a coupled nonlinear set of integro-differential equations connecting the single-particle GF to the self-energy 
operator is rather nontrivial. A characteristic feature of these equations is that, besides the single-particle GF, they 
involve also higher-order GF. The irreducible counterparts of the GFs, vertex functions, serve to identify correctly 
the self-energy as 

M^G^^ -G-^ 

The integral form of Dyson equation (|18|l gives A/ the physical meaning of a nonlocal and energy-dependent effective 
single-particle potential. This meaning can be verified for the exact self-energy through the diagrammatic expansion 
for the causal GF. 

It is important to note that for the retarded and advanced GFs, the notion of the proper part M — {P)p is symbolic 
in nature^. In a certain sense, it is possible to say that it is defined here by analogy with the irreducible many-particle 
T-matrix. Furthermore, by analogy with the diagrammatic technique, we can also introduce the proper part defined 
as a solution to the integral equation (|17|l . These analogues allow us to understand better the formal structure of 
the Dyson equation for the double-time thermal GF but only in a symbolic form . However, because of the identical 
form of the equations for GFs for all three types ( advanced, retarded, and causal ), we can convert in each stage of 
calculations to causal GF and, thereby, confirm the substantiated nature of definition H17|l ! We therefore should speak 
of an analogy of the Dyson equation. Hereafter, we drop this stipulating, since it does not cause any misunderstanding. 
In a sense, the IGF method is a variant of the Gram-Schmidt orthogonalization procedure. 

It should be emphasized that the scheme presented above gives just a general idea of the IGF method. A more exact 
explanation why one should not introduce the approximation already in P, instead of having to work out M , is given 
below when working out the application of the method to specific problems. 

The general philosophy of the IGF method is in the separation and identification of elastic scattering effects and 
inelastic ones. This latter point is quite often underestimated, and both effects are mixed. However, as far as the 
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right definition of quasi-particle damping is concerned, the separation of elastic and inelastic scattering processes is 
believed to be crucially important for many-body systems with complicated spectra and strong interaction. 
From a technical point of view, the elastic GMF renormalizations can exhibit quite a nontrivial structure. To obtain 
this structure correctly, one should construct the full GF from the complete algebra of relevant operators and develop 
a special projection procedure for higher-order GFs, in accordance with a given algebra. Then a natural question 
arises how to select the relevant set of operators {Ai,A2, ■■■An} describing the "relevant degrees of freedom". The 
above consideration suggests an intuitive and heuristic way to the suitable procedure as arising from an infinite chain 
of equations of motion H10|l . Let us consider the column 

A2 

where 



Ai^A, A2^[A,H], ^3 = [[Aff],ff],...^n = 



Then the most general possible Green function can be expressed as a matrix 

■A, 



\An) 



. ^1 ^2 



Ai)» 



This generalized Green function describes the one-, two- and n-particle dynamics. The equation of motion for it 
includes, as a particular case, the Dyson equation for single-particle Green function, and the Bethe-Salpeter equation 
which is the equation of motion for the two-particle Green function and which is an analogue of the Dyson equation, 
etc . The corresponding reduced equations should be extracted from the equation of motion for the generalized GF 
with the aid of the special techniques such as the projection method and similar techniques. This must be a final goal 
towards a real understanding of the true many-body dynamics. At this point, it is worthwhile to underline that the 
above discussion is a heuristic scheme only but not a straightforward recipe. The specific method of introducing the 
IGFs depends on the form of operators An , the type of the Hamiltonian, and conditions of the problem. 
Here a sketchy form of the IGF method is presented. The aim is to introduce the general scheme and to lay the 
groundwork for generalizations. We demonstrated in^ that the IGF method is a powerful tool for describing the 
quasi-particle excitation spectra, allowing a deeper understanding of elastic and inelastic quasi-particle scattering 
effects and the corresponding aspects of damping and finite lifetimes. In the present context, it provides an efficient 
tool for analysis of the mean fields and generalized mean fields of the complicated many-body models. 



IV. QUASI-PARTICLE DYNAMICS OF THE (sp - d) MODEL 



To describe self-consistently the spin dynamics of the extended sp — d model, one should take into account the 
full algebra of relevant operators of the suitable "spin modes" which are appropriate when the goal is to describe 
self-consistently the quasi-particle spectra of two interacting subsystem. 
We have two kinds of spin variables 



Let us consider the equations of motion 

[St,H,^a]- = -IN-' Y.^2S^-A^^P+il - ^^9(«pT«p+?T " 4i^p+<ii)] (20) 

pq 
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[5'„fc, ffs-d]- = -IN 1 Y^[2S^_gal^ap+g^ - Sf^_^{al^ap+q-'f - aJ^ap+^J] (21) 

pq 

Hs-d]- = -IN-' 5](5+_,4iap+,T - 5^_^4^ap+, J (22) 

pq 

[S+, Hd]- = Jq{SlSt^ ~ 5,^_,5+) (23) 

[5Z„i/,]_ =iV-V2^ „5.5-^^^^) (24) 

[aj^a^+fci, (£((? + fc) - e((7))aj^ag+fex + 

UN-'Y{al^al^p,^ap^aq+p,+ki - flg+p'T^-P'i^Pi^^+'^i) (^5) 
pp' 

[aj^aq+fei, = IN '/'^y^\Stp,{al^ap+p'^5p^q+k ~ a^piO'q+ki^q.p+p') 

pp' 

-Slp'{a\^o,p+p' i5p,q+k + a}p^aq+kiSq.p+p')] (26) 

From Ea. H20|) - Ea. l|26|l it follows that the localized and itinerant spin variables are coupled. Suitable algebra of 
relevant operators should be described by the 'spinor' (5') ("relevant degrees of freedom"), according to the IGF 
strategy. In principle, the complete algebra of the relevant "spin modes" should include the longitudinal components 
cr^ and 51 . However, the correlation of the longitudinal spin components are rather small at low temperatures and 
become essential with the approaching to the Curie temperature. The calculation of the Green function for the 
longitudinal spin components is a special non-trivial task^. Since we are interesting here in the low-energy spin- wave 
type of excitations, we will consider the transversal components only. 

The model Hamiltonian H — Hg+Hs-d + Hd was used in^;^ for calculations of the spin- wave spectra and was called 
the modified Zener model. In this model, as applied to transition metals, the itinerant electrons are described by a 
Hubbard Hamiltonian and the itinerant electron couples the localized spin (Hund's rule coupling) by a term -ffs-d- 
Because of the inequivalent spin systems, localized and itinerant, a consequence of the model is the existence of acoustic 
and optic branches of the quasi-particle spectrum of spin excitations. In DMS the local antiferromagnetic interaction 
Hs-d produces the coupling between the carriers ( which are holes in GaMnAs ) and the Mn magnetic moments ( s 
— 5/2) which leads to ferromagnetic ordering of Mn spins in a certain range of concentration. The Kondo physics is 
irrelevant in this case, but the fully determined and consistent microscopic mechanism of the ferromagnetic ordering 
is still under debates^i*^^. An important question in this context is the self-consistent picture of the quasi-particle 
many-body dynamics which takes into account the complex structure of the spectra. 

A. Spin Dynamics of the s — d Model. Scattering Regime. 

In this section, we discuss the spectrum of spin excitations in the sp ~ d model. We consider the double-time 
thermal GF of localized spins'*^ which is defined as 

G+-{k; t - t') =« S+{t),SZk{t') »= -iO(t - t') < [S+{t), SZkit')]- >= 

/ + 00 
dujexp{~iLut)G+-{k;uj) (27) 
-oo 

The next step is to write down the equation of motion for the GF. 

Our attention will be focused on spin dynamics of the model. To describe self-consistently of the spin dynamics of 
the sp^d model, one should take into account the full algebra of relevant operators of the suitable "spin modes" which 
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are appropriate when the goal is to describe self-consistently the quasi-particle spectra of two interacting subsystems. 
We introduce the generahzed matrix GF of the form 



= G{k;uj) (28) 



Here 

Equivalently, we can do the calculations with the matrix of the form 

{ <<^^-l^-» «sM.f.^» (29) 

but the form of Ea. (|28|l is slightly more convenient. 

Let us consider the equation of motion for the GF G{k;uj). By differentiation of the GF << {t)\B{t') >> with 
respect to the first time, t, we find 

.«5J|B (30, 

pq 

+iV-i/2 ^ « (5|5+_^ - 5,%5+)|i? (31) 

where 

B 



'-k 



Let us introduce by definition irreducible {ir) operators as 



(S"!) — S*^- < S'q > 5q^Q] Wp+qa^pa) — Op+qcr^pCT— < 0.\a^pa > <^g,0 (32) 



((5|) - (S^k-q) = ((5|) - (5,^_,) ^'•5+) - (0, - ^fc_,)5+ (33) 



-1 

From the condition (|12|l 

< [((5|) "Sl_^ - {Sl_^) -5+ - (^, - ^fc-,)5+), 51,] _ >= 

one can find 

2 < S-o^ > 



(34) 



^<[S\)^^{Sl)^^ >■ K-+^<SZqS+> (35) 

Using the definition of the irreducible parts the equation of motion Ea. pOl ) can be exactly transformed to the following 
form: 

r!i « S+\B »^ +^2 « (7+\B »^= |^~^^"'|+ « MB »^ (36) 



where 



fii = - ^^(Jo - Jk) - N-''^Y.^Jq J^_,)^2__^ _ /(„^ _ (37) 
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]Sf ^ ~' "I'^-'i" ' TV 

9 9 



= Y^{n^+ni); 0<n<2 



n 



The many-particle operator Ai reads 

pg 

+iV-V2 ^ - (SLJ '''5+) (39) 

and it satisfies the conditions 

< [Ai,SZq]- >=< [Ax,aZq]- >=0 

To write down the equation of motion for the Fourier transform of the GF << ct^ (t), B{t') », we need the auxiliary 
equation of motion for the GF of the form 

« al^ap+ki{t),B{t') » 

For this we have to write the equation of motion for it after differentiation with respect to the first time variable 
t and extract the corresponding irreducible parts. Then, we obtain, after the Fourier transformation, the following 
equation: 

(w + e{p) - e{p + k)- 2/7V-V2 < 5^ > -U{n^ - n^)) « a^^ttp+^jB + (40) 
UN-\fp^ - fp+ki) « (T+\B +IN-'/\fp^ - fp+ki) « S+\B »„= 





qr 

UN-^ X « {al^a\^^^aq^ap+r+ki - ap+rT^g-ri^giOp+fci) 'H-^ 



qr 



Let us use the following notation: 

A2 = -IN~'^/'^'^[Sl^{al-faq+r^6p+k,q - al^ap+ki6p,q+r) - 

qr 

(Sl^) '''{al^ap+kiSp,q+r + al^aq+ridp+k,q)] + (41) 

qr 

uJp,k = (w + e{p) - e{p + k)-A) (42) 
A = 277V-V2 <S^> -C/(n| - n^) = 2IS -Um = Aj + Au (43) 

ro{k,u;) = iV^i J2 ^^'^''"^"'^ (44) 
p '^p,k 
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Now we consider the GF << (T^(t),B{t') » . Similarly to Eq.(|2ni), we have 

- N^/^Ixl{k,Lo) « S+\B »^ +(1 - UxUk^Lo)) « a+\B »^= 



mIi. ^)+E— (45) 

Here the following definition of the irreducible part for the Coulomb correlation term was used 

(t 1' 1^ t ^ 

t X t 

~ < "g-ri^-ji > ^q-rWp+r^(^P+kl (46) 

The operator A2 satisfies the conditions 

< [A2,SZk]- >-< [A2,'JZk]- >=0 
In the matrix notation the full equation of motion for the GF G{k:Lu) can now be summarized in the following form: 

ClG{k;to) ^ i + Y,Hp)Dip;to) (47) 
p 

(G(fc;^))t(17)t = (/)t + ^(^(p;c.))t($(p))t 



where 



Dip;u;)={ / *(P)= n "M (49) 

To calculate the higher order GFs in f l47|) . we differentiate its r.h.s. with respect to the second-time variable (f). Let 
us give explicitly one of the four equations. After introducing the irreducible parts as discussed above we get 



<<Ai|S'.7>> << ^ilfTl.. >>\ , 



N 
p'q' 



Here the symbolic notation for the three equation of motions were used with « — 1, 2,3. The quantity Ai in the l.h.s. 
of (EH) should be substituted by 



A, = ((5^) -^+_^ - (S^^^^) "-3+) - 
A3 = 2S'|_^a],^ap+,x 

In the matrix notation the full equation of motion for the GF D{k;uj) can now be written in the following form: 

i77)(p;a;)=^<I.(p')I)i(p';c^) (51) 
p' 



where 



«A2\A\» «AMl»' ^ ' 
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Combining both (the first- and second-time differentiated) equations of motion, we get the " exact" ( no approximation 
has been made till now) "scattering" equation 

nG{k;io) = i + J2Hp)PiP,PMp'm^)~' (53) 
pp' 

This equation can be identically transformed to the standard form Eq. ( I16|l 

G = Go + GoiY. r^<^{p)P{p,p'Mp')i-^)Go 
pp' 

G = Go + GoPGo (54) 
Here we have introduced the generalized mean- field (GMF) GF Go, according to the following definition: 

Go = ^'^I (55) 

The scattering operator P has the form 

P = /-i^<|(p)P(p,p')|.(p')/-i (56) 
pp' 

Here we have used the obvious notation 

U<A2|At » «A,\Al») ^^^^ 
As is shown above, Ea.(' l54|l can be transformed exactly into the Dyson equation f I18|l 

G = Go + GoA/Go (58) 

with the self-energy operator M given as 

M = {P)P (59) 
Hence, the determination of the full GF G has been reduced to that of Go and M. 

V. GENERALIZED MEAN-FIELD GF 

From the definition ( I55|) , the GF matrix in the generalized mean-field approximation reads 

where 

i? = (1 - ux^)ni + n^iN^'^i 

Let us write down explicitly the diagonal matrix elements G}^ and G^ 

9 ^ 

where 

X^k,u;) = Xoik,^){l - Uxlik^u))-^ (63) 
S = N-^'^ < > 

To clarify the functional structure of the generalized mean-field GFs (|61l ) and Ht)2l ) , let us consider a few limiting 
cases. 
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A. Uncoupled Subsystems 

To clarify the calculation of quasiparticle spectra of coupled localized and itinerant subsystems, it is instructive to 
consider an artificial limit of uncoupled subsystems. We then assume that the local exchange parameter 1 = 0. In 
this limiting case we have 

25 

« Sk \s_, » = ^ _ -^j^ _ _ _^ ^^^j^ _ j^_^)(2i^- + k;[+) ^^^^ 

«a+\aZ,»'=rik,u;) (65) 

The spectrum of quasi-particle excitations of localized spins without damping follows from the poles of the generalized 
mean-field GF ^ 

u{k) = 5(Jo ~Jk) + ~ Y.(J, - Vfc)(2i^r + K^^) (66) 

q 

It is seen that due to the correct definition of generalized mean fields we get the result for the localized spin Heisenberg 
subsystem which includes both the simplest spin-wave result and the result of Tyablikov decoupling as limiting cases. 
In the hydrodynamic limit fc — *■ 0, a; — > it leads to the dispersion law Lo{k) — Dk^. 
The exchange integral Jk can be written in the following way: 

Jfc =^exp(-ifci?,)J(|i?,|) (67) 

i 

The expansion in small k gives 



i i i 

Here r? = k/k is the unit vector. The values Jk-q can be evaluated in a similar way 

Jk-g ^Jg- {kVg)jg + ^{kVgfjg + ■■■ (69) 

ikVg)Jg = ~i"^{kR,)J{\R,\)exp{~iqR,) 

i 

{kVgfjg = -]^Y.^kRifj{\R,\)e^^{-iqRi) 

i 

Combining Eq.lEHIl, Eq.®, and Eq.ljnHI) we get 

« S+\SZk »°= ^ (70) 

uj-uj{k) 

^(k ^ 0) = (5( Jo - Jk) + ~ Y.(Jg ~ Jg-k){2Kl^ + K-+)) ^ 



2NS ■ 
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(1^0 + j5E^.(2^r+^r 

V-g = ^(fci?^)V(|i?,|) exp (-i(fi?. 



Let us now consider the spin susceptibility of itinerant carriers Eq. (|65h in the hydrodynamic limit fc ^ 0, a; ^ 0. It 
is convenient to consider the static limit of Eq. (|65l ) 



<< '^k r-k » I'^^o - TT.^sn. n\ i'l) 



xg(fc,o) 

1-C/Xg(fc,0) 

Xoik, 0) = i E ,(g) ^^'e{q + k)- Au 
Ajj = U{ni — ni) = Um 
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To proceed, we make a small-k expansion of the form 



e(g + k) -e{q) = (fcV,)e(g) + i {k'^,)h{q) + ■■■ (72) 



U q ^ q 

The poles of the spin susceptibihty of itinerant carriers are determined by the equation 

l-Uxl{k,Lo)^Q (73) 

In another form this reads in detail 

^ ST^ ^ fq+kl 



N ^ e{k + q)- e{q) + Ay - w 
If we set Lo = E{k) and put then fc = 0, we get the equation for the excitation energy E{k = 0) 



U V — > 

1= /vE 



N^Au-E{k^O) Au-E{k^O)U 

which is satisfied if E{k = 0) = 0. Thus, a solution of Ea. l|73l ) exists which has the property Wmk^i^ E{k) = and 
this solution corresponds to an acoustic spin-wave branch of excitations 

E{k) = - E(/?T + .Ui){kVqfe{q) + ^ ^(/^T ^ hi){ky A^l)? (74) 

^ q U q 

to = eik + q) -e{q) + Au 

It is seen that the stiffness constant D2 can be interpreted as expanded in For the tight-binding electrons in s.c. 
lattice the spin wave dispersion relation D2k^ becomes 

= (3(n, -n,))-i5][i^^-Ai)|v,e(g)p - ihl±Myl,^^] = (75) 

q ^ 

(3(n| - "i))"^(— r — E(At - fil){kx sm{qxa) + kySm{qya) + k^ siniq^a))'^ - 
Au Y 



ta'^ ^ifil + fql ) i^l cos q^a + ky cos qyU + k^ cos q^a)) 



B. Coupled Subsystems 

The next stage in the analysis of the quasi-particle spectra of the (sp — d) model is the introduction of the nonzero 
coupling I. The full generalized mean field GFs can be rewritten as 

2 9 

<< >>°= = - = (76) 

' - Im- S{J„ - Jk) - ^EqiJ, - Jq-k){2K^^' + Kq+) +2PSr{k,iu) 



Here the notation is used 



2PS 

Ueff ^ ^ ~ ^ _ ' ^=(^^-^1) 

The expression Eq.|(77|) coincides with that for the itinerant spin susceptibility as calculated in^. It is instructive to 
consider separately the four different cases, 
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(i) / 7^ 0, J = 0, t/ = 

(ii) I ^0,J ^0,U = 

(iii) I ^0,J = 0M ^0 

(iv) I ^0,J ^0,U ^0 



The first case I 0, J = Q,U = corresponds to a model which is commonly called the Kondo lattice model. It 
can be seen that GFs (|7(jl ) and (|77l ) are then equal to 

9 S' 

« S+\S-^ „ (78) 



In order to calculate the acoustic pole of the GF {TSl), we make use of the small (fc, w) expansion. Hence we get 

«S+\S-_k»'^ (80) 

2S' 



25(l + ||)-i 



It follows from Eg. 1181)1 ) that the stiffness constant D is proportional to the total magnetization of the system. 

2. 

In the second case / ^ 0, J 7^ 0, ?7 = 0, we get 

g+\g^ >>"= — (81) 

' Lo-Im-S{.h-Jk)-^Y.q{Jq-J,-k){2K^q'-+K^+) + 'iPSxl{k,oj) 

<< CTfc |o-_fc >> 2j2g — : (82) 

1-I7^4(fc,^) 

In order to calculate the acoustic pole of the GF H81I '). we make use of the small (fc, oj) expansion again. We then get 

« 5+1^1, »"« (83) 

2^(1 + ^)-! 

^ - (1 + %)-^D,B - (1 + t)"M2FAf E,(/,T + U){k^,Ye{q) - ^ E,(/,T " /gi)(fcV,6(gO)2] 

It follows from Eas. (|80l 1 and (|83l ) that the stiffness constant D is proportional to the total magnetization of the 
system. 



The third case / t^O, J = 0, corresponds to a model which is called the modified Zener lattice model^. It 

can be seen that in this case GFs (1761 ) and (|77l ) are equal to 



2S_ 

uj - Im + 2PSx''{k,uj) 



« S+\SZ, »"= . , ... (84) 
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The results obtained here coincide with those of Bartel^^. The excitation energies for the locahzed spin and spin 
densities of itinerant carriers are found from the zeros of the denominators of << S^\SZk >>° and << (j'^\aZk >>" 
which yield identical excitation spectra, consisting of three branches, the acoustic spin wave i?'"^(fc), the optical spin 
wave E°P{k), and the Stoner continuum E^^{k) 

E'"'{k) = Dk^ 



JJE°P 

E"P{k) = E°P - Dil - -j^)k'; E°^ = /(m + 2^) 



E^\k) ^ e{k + q) ^ e{q) + A 



4- 

The most general is the forth case, I ^ 0, J Q,U ^0. The total GF of the coupled system is given by Eq.|7Bl). 
The magnetic excitation spectrum follows from the poles of the GF ( I6U|I 

R={1- Ux'o)^i + n^iN^'^i = 

and consists of three branches - the acoustic spin wave E°''^{k), the optical spin wave E°P{k), and the Stoner contimmm 
E'^^k). 

Let us, as a first approximation, consider the last term in its denominator , which is the dynamic spin susceptibility 
of itinerant carriers, in the static limit, without any frequency dependence. The GF Ea. (|76h then becomes equal to 



+ 25 

« Sk \S_k » ~ _ j^) _ _1_ ^^(j^ _ Jg-fe)(2X- + Kg+) + 2PSr{k,0) ^^^^ 

It is possible to verify that in the limit fc ^ 

2l'Sx%K 0) « - ^ E(/«T + U)ik^, + E^Z-^T - Ui){kV,e{q)r (87) 

9 9 

Then for w, fc ^ Eq.(|SH|) becomes 

«S+\SZk»"^ (88) 

2S 

u; - D,k^ - ^ E,(/9T + /,i)(fcV,)2e(g) + ^ Y.,{U^ ~ UKkyAo))' 

This expression can be expected to be qualitatively correct in spite of the primitive approximation. The spectrum of 
Stoner excitations is given by 

E^\k) ^ e{k + q) - eiq) + A (89) 

In addition to the acoustic branch there is an optical branch of spin excitations. This can be seen from the following: 
For fc = we get for R = the quadratic equation in u with two solutions, lu = and u = I(rn + 2S) = E^^ . In the 
hydrodynamic limit, k ^ 0, uj ^ the GF H75|l can be written as 



25" 

u -~E^'=(fc) 

where the acoustic spin wave energies are given by 



«S+\SZ,»"^—^^ (90) 



£;-(fc) = Dk' = (|[^o + ^ E^'/(2i^f + K^+)] (91) 



4:^E(/'^T + Min^,rei<l) + T^;^E(/'^T - /,i)(r?V,6(g))%2 



2iV25'^'^' ^' ^ ' NA2S 

1 q 
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For the optical spin wave branch the estimations can be carried out as in paper— 

E°P{k) = E°P ~ D°Pk^ (92) 

In the GMF approximation the density of itinerant electrons ( and the band splitting A) can be evaluated by solving 
the equation 

= ^ ^[exp(/3(e(fc) + -IS- ep)) + 1]-' (93) 

k 

Hence, the stiffness constant D can be expressed by the parameters of the sp — d model Hamiltonian . 

VI. EFFECTS OF DISORDER IN DMS 

We now proceed to give a simple and qualitative discussion of the effects of disorder in DMS to give just a flawor 
of ideas how the disorder can be included in the IGF scheme. The full threatment of disorder effects require the 
consideration of damping effects and will be considered separately. 

The main aim of the investigation of DMS is to give a successful microscopic picture of the ferromagnetic ordering of 
localized spins induced by the interaction with the spin density of itinerant charge carriers. As it was stated above, 
a suitable model, which may be used for investigation of this problem ( at least at the initial stage ) is a modified 
Kondo lattice model (|H1) 

ij a i 

Here Vi projects out sites occupied by Mn atoms, i.e.: 

_ J 1 if site i is occupied by Mn 
* [ if site i is occupied by Ga 

This model is relevant for the doped II- VI or III-V compound. The essential feature of the model is that it describes 
a mechanism of how the spins of carriers (electrons or holes ) become polarized due to the local antiferromagnetic 
exchange interactions with localized spins. In A^-^Ji^MrixB^ the main magnetic interaction is an antiferromagnetic 
exchange between the Mn spins and the charge-carrier spins. The superexchange term — JijSiSj is 

antiferromagnetic also but is as a rule rather small in the concentration range of interest ( x ~ 0.05). In the case of 
Mn-doped III-V compounds the antiferromagnetic superexchange interaction will generally reduce the ferromagnetic 
ordering temperature. As a result, the carrier-induced ferromagnetism in DMS arises due to the effective ferromagnetic 
interaction between the Mn spins. In other words, the ferromagnetism in this system is most probably related to the 
uncompensated Mn spins and is mediated by holes. The density of Mn ions cmu is greater than the hole density 
P, CMn S> p. The optimal interrelation of both the magnitudes is a delicate and subtle question and was analyzed 
recently in paper—. It was shown that the concentration of free holes and ferromagnetically active Mn spins were 
governed by the position of the Fermi level which controls the formation energy of compensating interstitial Mn 
donors. The experimental evidence has been provided that the upper limit of the Curie temperature is caused by 
Fermi-level-induced hole saturation. In order to provide a suitable treatment of the spin quasiparticle dynamics it 
is necessary to take into account the effects of disorder since the Mn ions are assumed to be distributed randomly 
with concentration c. This is positional disorder. There is variation of site-energy of nonmagnetic origin due to 
the substitution of A atom with Mn ion. The detailed nature of the disorder is not fully clear. In paper^Si, it was 
shown that the dominant fraction of the Mn atoms are on either substitutional sites or specific sites shadowed by 
the host atoms. This reveals that the majority of the Mn atoms are on specific ( nonrandom ) sites commensurate 
with the lattice, but that does not necessarily imply that all of the Mn atoms are in substitutional positions. For 
X > 0.05 an increasing fraction of Mn spins do not participate in ferromagnetism. It can be related with an increase 
in the concentration of Mn interstitials accompanied by a reduction of Tc- There are indications of an increase in 
Mn atoms in the form of random clusters not commensurate with the GaAs lattice. However, these results require 
the independent confirmation. The conclusion that there is a maximum in Tc due to the Fermi level pinning is a 
conjection only. There are evidences that the largest values of Tc have been found to be considerably larger than 110 

■|^^21j22j35 

It follows from Eq. (|94l ) that the spin dynamics of a modified KLM will be described by the GFs in the lattice site 
representation for a given configuration 

« S+\SJ- » « a+\aj » 
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and instead of Eq. (|28|l the lattice GF should be considered 

''^'i^^^^ «t^">=^UG.,M (95) 
«a+\S^ » «c7+\a^ »J ^ ^ ^ 

In order to provide a simultaneous and self-consistent treatment of the quasi-particle dynamics including the effects 
of disorder, a sophisticated description of disorder should be done. Most treatments remove disorder by making a 
virtual-crystal-like approximation in which the Mn ion distribution is replaced by a continuum. A more sophisticated 
approach for treating the positional disorder of the magnetic impurities inside the host semiconductor is the CPA — . 
The CPA replaces the initial Hamiltonian of disordered system by an effective one which is assumed to produce no 
further scattering^. It describes reasonably well the state of itinerant charge scattering in disordered substitutional 
alloys Ai^.j;Bx. 

In order to simplify the discusson here, we will deal with a much simpler and less sophisticated description. The 
approximation discussed below should be considered as a first, crude approximation to a theory of disorder effects 
in DMS. Since the detailed nature of disorder in DMS is not yet established completely, we will confine ourselves 
to the simplest possible approximation. Let us remind that the IGF method is based on the suitable definition of 
the generalized mean fields—. To demonstrate the flexibility of the IGF method, we show below how the mean field 
should be redefined to include the disorder in an effective way. The previous definition of the irreducible spin operator, 
Eq.(|321), should be replaced by 

{Sg) = Sq - c< >Sqfi; (a^+^^flp^) " = a^+^^ap^- < al^ap„ > Sq^ (96) 

Here < Sz >= N~^/^Sz corresponds to the configuration average. The average < Sz > denotes the mean value of 
for a given configuration of all the spins. We omitted here the variation of site energy of nonmagnetic origin. 
The consequences of this choice manifest themselves. It means precisely that in a random system the mean field is 
weaker as compared to a regular system. The approximation is conceptually as simple as an ordinary mean field 
approximation and corresponds to the virtual crystal approximation. The situation is then completely analogous to 
the previous one considered in the preceding sections. We get for the configurationally averaged GFs 

2cSz 

Lo - Im + 2PcSzXo{k,oj) 



« St\Sr »o ^« S+\SZ, »°« , (97) 



« at\aj »o =« at\aZ, »"« .l^gj'l. . (98) 



j-Ir 



These simple results are fully tractable and are the reason for their derivation. 

It is worth to note that in the case of the modified Zener model which contains the correlation (Hubbard) term, the 
effects of disorder should be considered on the basis of a similar model^ 

i7 = ^ ^ Ujal^a-ja + lyiU^^n^i - ^ 2IviS^Si (99) 

ij (7 i i 

The Coulomb repulsion is assumed to exist only on lattice sites occupied at random by Mn atoms. The approach mostly 
usedS to calculate stiffness constant within a random version of the Hubbard model was based on the random phase 
approximation, where the electron-electron approximation was taken into account in the Hartree-Fock approximation 
and the disorder in the CPA. It is therefore very probable that within this approach the formation of magnetic clusters 
can be reproduced; the formation of the clusters is thus strongly enviromcntal-dcpcndcnt. However, the calculation 
of the spatial GF Ea.f l95|) . for the model, Ea.f l99|) . is rather a long and nontrivial task and we must avoid considering 
this problem here. We hope, nevertheless, that the description of the disorder effects, as given above, gives a good 
first approximation as far as the the irreducible Green functions method is concerned. A more detail consideration of 
the state of itinerant carriers in DMS including a more sophisticated treatment of disorder effects will be carried out 
separately. 



VII. CONCLUSIONS 



In summary, we have presented an analytical approach for treating the spin quasi-particle dynamics of the gener- 
alized spin-fermion model, which provides a basis for description of the physical properties of magnetic and diluted 
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magnetic semiconductors. We have investigated the influence of the correlation and exchange effects for interact- 
ing systems of itinerant carriers and locahzed spins using the ideas of quantum field theory for interacting electron 
and spin systems on a lattice. The workable and self-consistent IGF approach to the decoupling problem for the 
equation-of-motion method for double-time temperature Green functions has been presented. The main achievement 
of this formulation was the derivation of the Dyson equation for double-time retarded Green functions instead of 
causal ones. That formulation permits to unify convenient analytical properties of retarded and advanced GF and 
the formal solution of the Dyson equation which, in spite of the required approximations for the self-energy, provides 
the correct functional structure of single-particle GF. The main advantage of the mathematical formalism is brought 
out by showing how elastic scattering corrections (generalized mean fields) and inelastic scattering effects (damping 
and finite lifetimes) could be self-consistently incorporated in a general and compact manner. In this paper, we have 
confined ourselves to the elastic scattering corrections and have not considered the damping effects. This approach 
gives a workable scheme for definition of relevant generalized mean fields written in terms of appropriate correlators. 
A comparative study of real many-body dynamics of the generalized spin-fermion model is important to characterize 
the true quasi-particle excitations and the role of magnetic correlations. It was shown that the magnetic dynamics 
of the generalized spin-fermion model can be understood in terms of combined dynamics of itinerant carriers, and of 
localized spins and magnetic correlations of various nature. The two other principal distinctive features of our calcu- 
lation were first, the use of correct analytic definition of the relevant generalized mean fields, and second, the explicit 
calculation of the spin-wave quasiparticle spectra and its analysis for the two interacting subsystems. This analysis 
includes all of the interaction terms that can contribute to the essential physics. Thus the present consideration is the 
most complete analysis of the quasiparticle spectra of the spin-fermion model of magnetism within the generalized 
mean field approximation. These applications illustrate some of subtle details of the IGF approach and exhibit their 
physical significance in a representative form. 

As it is seen, this treatment has advantages in comparison with the standard methods of decoupling of higher order 
GFs within the equation-of-motion approach, namely, the following: At the mean-field level, the GF one obtains, is 
richer than that following from the standard procedures. The generalized mean fields represent all elastic scattering 
renormalizations in a compact form. 

The approximations ( the decoupling ) are introduced at a later stage with respect to other methods, i.e., only into 
the rigorously obtained self-energy. 

The physical picture of elastic and inelastic scattering processes in the interacting many-particle systems is clearly 
seen at every stage of calculations, which is not the case with the standard methods of decoupling. 
Many results of the previous works are reproduced mathematically more simply. 

The main advantage of the whole method is the possibility of a self- consistent description of quasi-particle spectra and 
their damping in a unified and coherent fashion. However, in the present paper, for the sake of clarity, we concentrated 
on the clear presentation of the quasiparticle many-body dynamics within a generalized mean field approximation. 
This explains why we confine ourselves by consideration of disorder effects in the simplest VGA. The consideration 
of disorder effects beyond VGA includes many intrinsic specific problems and deserves a separate investigation. The 
irreducible GFs methods will be generalized to treat these problems in separate publications. 

Thus, this new picture of an interacting spin-fermion system on a lattice is far richer and gives more possibilities for 
analysis of phenomena which can actually take place. In this sense, the approach we suggest produces more advanced 
physical picture of the quasi-particle many-body dynamics. Our main results reveal the fundamental importance 
of the adequate definition of generalized mean fields at finite temperatures, which results in a deeper insight into 
the nature of quasi-particle states of the correlated lattice fermions and spins. The key to an understanding of the 
situation in DMS lies in the right description of the interplay of interactions and disorder effects for coupled spin and 
charge subsystems. Consequently, it is crucial that the correct functional structure of generalized mean fields was 
calculated in a closed and compact form. The detailed consideration of the state of itinerant charge carriers in DMS 
along this line will be considered separately. 
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